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Abstract 

The Quadrantid meteor shower is among the strongest annual meteor showers, and has drawn the attention of scientists for several decades. The stream is unusual, 
among others, for several reaso ns: its very shor t duration around maximum activity (« 12 - 14 hours) as detected by visual, photographic and radar observations, its 
recent onset (around 1835 AD jOueteletl . ll839[) ) and because it had been the only major stream without an obvious pai'ent body until 2003. Ever since, there have 
been debates as to the age of the stream and the nature of its proposed parent body, asteroid 2003 EHi. 

In this work, we present results on the most probable age and formation mechanism of the narrow portion of the Quadrantid meteoroid stream. For the first time 
we use data on eight high precision photographic Quadrantids, equivalent to gram - kilogram size, to constrain the most likely age of the core of the stream. Out of 
eight high-precision photographic Quadrantids, five pertain directly to the naiTow portion of the stream. In addition, we also use data on five high-precision radar 
Quadrantids, observed within the peak of the shower. 

We performed backwards numerical integrations of the equations of motion of a large number of ’clones’ of both, the eight high-precision photographic and 
five radai* Quadrantid meteors, along with the proposed parent body, 2003 EH]. According to our results, from the backward integrations, the most likely age of 
the narrow structure of the Quadrantids is between 200 - 300 years. These presumed ejection epochs, corresponding to 1700 - 1800 AD, are then used for forwai'd 
integrations of large numbers of hypothetical meteoroids, ejected from the parent 2003 EHi, until the present epoch. The aim is to constrain whether the core of the 
Quadrantid meteoroid stream is consistent with a previously proposed relatively young age (% 200 years). 
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1. Introduction 
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The Quadrantids are among the most active meteor showers, 
^ reaching a peak activity of Zenithal Hou rly Rate (ZHR) ~ 110 
—I - 130 on 3-4 Januar y each year ( Shelton . 19651 Hindlevl 1970 : 

^ iHughes and Tavlor . 1977h . as determined by photographic, vi¬ 
sual, video and radar technique s. The stream has recently been 
linked with asteroid 2003 EHi ( Jenniskeni, 2004). 

The Quadrantid shower is unusual among meteoroid streams 
presently visible at the Earth for several reasons. Eirstly, the 
Quadrantid meteor shower has a short duration of maximum ac¬ 
tivity, which we will call the ’core’ or the ’’narrow structure” of 
the stream. The Eull-Width at Half-Maxir num (lAVHM) of the 
core activity is ^ 0.6 days ( Shelton . 19651 Hughes and Tavloi 
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Il977 : Browei . 20061) for visual-sized particles, implying that 
this central portion is very young, while the shower as a whole 
has an overall duration of significant length ~ 4 days. Sec¬ 
ondly, iyias_onlybecomeactiyej^ently, being first noted circa 
1835 ( Quetelel , 18391 Eishei , 1930l) . Moreover, the activity of 
the shower has changed dramatically over the last 150 years, 
from a very weak s hower to among the strongest visible at the 
Earth djenniskensl 20061) . Einally, recent radar observations 
( Brown et al. , l2010bl) suggest low level activity of the shower 
persisting for a few months (November to mid January), sug¬ 
gesting the stream has an older component as well. 

Presently, the presumed parent body of the core of the Quad¬ 
rantids is the Near Earth Object (NEO) 2003 EHi. 6 Through¬ 
out this paper we will refer to it as asteroid 2003 EHi. The 
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object has been classified as an Amor type asteroid, although 
its nature is arguable based on dynamical criteria. Asteroid 
2003 EHi has a short-period comet-like orbit, with a Tisserand 
parameter with respect to Jupiter of Ty « 2.0, but currently 
shows no evidence of cometary activity, suggesting that it is a 
strong candidate fo r either a recently dormant or extinct comet 


dKoten et alll2006l) . 


Prior to the discovery of 2003 EHi, a few other ob¬ 
jects with less similar orbits had been connected to 
the Quadrantid meteoroid stre a m, m o st notably comet 


Babadzhanov and Obrubov 

1992) and comet C/1490 Y1 

dHasegawal 19791 iLee et al. 

20091 

Williams and Wu. 1993; 


ship of these bodies to the Quadrantids remains unclear. 

The earliest attempt to tackle th e age of the Quadrantid me 


teoroids stream can be attributed to Hamid and Youssefl ( 1963 ). 


The authors carried out a numerical secular perturbation anal¬ 
ysis on the orbit of six doubly photographed Quadrantids and 
discovered large variations in the perihelion distance and the in¬ 
clination of the stream orbit, with a period of 4000 years. Based 
on the backward secular solutions, the authors argued that the 
orbital elements of the six meteors were similar around 3000 
years ago. 

Williams et al. ( 19791) calculated the secular variations of 


the orbital elements of the mean Quadrantid stream and con¬ 
cluded that the Quadrantid meteoroid stream may have resulted 
from two major comet break-ups about 1690 and 1300 years 
in the past, where the resulting meteoroids converged into their 
present orbit around 200 - 150 years ago, explaining the recent 
appearance of the stream. Similar work was also performed by 
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Hughes et al.l (ll97^. 

Hasegawa ( 19791) . was the first to propose a potential parent 


body for the Quadrantids, noticing a similarity between the or¬ 
bits of the mean Quadrantid stream and comet 14911 (= C/1490 
Yl), recorded in ancient Chinese observations. However, only 
a parabolic solution was assumed for the orbit of comet 1491 
I, due to the low observational accuracy in the position of the 
comet. Based on the arguable similarity between the orbits of 
1491 I and the Quadrantids, the authors concluded that 1491 I 
had been a periodic Jupiter-family comet, which suffered a very 
close encounter with Jupiter and was perturbed into a longer 
period orbit, where the orbital association with its meteoroid 
stream was lost. 

Assuming that 1491 I was the actual parent of the Quad- 
rantids using the calcul ated orbital elements of the comet. 


Williams and Wul (Il993h concluded that the stream was cre¬ 


ated ~ 5000 years ago. Based on the hypothesis of a very 
shallow close en c ounter between comet 1491 I and Jupiter, 


Williams and Wul (119931) demonstrated that prior to the en¬ 


counter with Jupiter, the eccentricity of the orbit of the comet 
must had been e 0.77. The newly derived value for the ec¬ 
centricity was used for backwards integration of the orbit of 
the comet to about 5000 years. Then, that epoch was used for 
the meteoroids ejection whose orbits were integrated forward 
until the present. The authors argued that the observed mean 
orbital elements of the stream is consistent with dust particle 
ejection ~ 5000 years ago. However, the lack of precise orbital 
elements for 14911, along with a hypothesized close encounter 
with Jupiter, rend ers the later conc lusion uncertain. For a simi¬ 
lar work, see also Lee et alj ( 2009 ). 

Another possible parent of the Qua drantids is the c omet 
96P/Machholz (formerly P/1986 VIII). iMcIntoshI (Il990l) cal¬ 
culated the secular precession of the orbital elements for the 
Quadrantids and comet 96P/Machholz and found that the long¬ 
term evolution of both orbits is strikingly similar, except for 
their precession cycles being shifted by a period of 2000 years. 
The author suggested that the stream was quite old and the 
phase shift in the precession cycles is due to the differential pre¬ 
cession of the orbits of the stream and the comet. Moreover, the 
author argued that Quadrantids may be a part of a larger com- 
pl ex of meteoroid stream, belonging t o comet 96P/Machholz. 


Babadzhanov and Obrubov (1 19921) integrated the orbits of 


three test particles similar to that of comet 96P/Machholz for 
8000 years back in time. Then 20 test particles were ejected 
from the nucleus of 96P at the epoch of 4500 B.C and inte¬ 
grated forward until 3000 AD. For that period of 7500 years, 
the authors argued that meteoroids released by 96P can pro¬ 
duce eight meteor showers on Earth within one precession cy¬ 
cle of the argument of perihelion w of the meteoroids. Six of 
these showers have been identified as; the Quadrantids, the Ur- 
sids. Southern 6 - Aquarids, daytime Arietids, Carinids and a - 
Cetids. That led the authors to conclude that 96P/Machholz is 
the parent of the Quadrantid mete oroid stream. For addit ional 
and more extensive work, see al so Bones and Jones! (Il993h and 
Kahucho va and N eslusanI (2007 ). 

Jenniskens et al.l ( 1997 ) used ~ 35 doubly-photographed 


Quadrantids taken in 1995 by the Dutch Meteor Society (DMS) 


(iBetlem et allll995() . to argue that the age of the central portion 
of the Quadrantid stream is only ~ 500 years old and the parent 
may be hidden on a Near Earth Obje ct (NEO) - like orbit. With 
the d i scovery of 2003 EH| in 2 003 (Jenniskens and Marsden . 
2 OO 3 I: McCluskv et al. . 2003h . Jenniskens (2004) noted the 
striking similarity between the current orbit of the Quadrantids 
and the orbit of 2003 EHi and proposed a sibling relationship. 


Wiegert and BrownI (120051) estimated an approximate age of 


200 years for the core of the Quadrantids, based on the nodal 
regression rate of the stream and forward integration of mete¬ 
oroids, released by 2003 EHi circa 1800 AD. The authors con¬ 
cluded that meteoroids released prior to 1800 AD appear on 
the sky at much earlier times than the first reported appearance 
around 1835. 


The main goal of this work is to estimate the most probable 
age of the central portion of the Quadrantid meteoroid stream 
and its mode of formation (e.g. cometary sublimation vs. as- 
teroidal disruption). We seek to first constrain the approximate 
formation age assuming 2003EH1 is the parent by first perform¬ 
ing backward integrations of high precision Quadrantid mete¬ 
oroids to compare the orbital similarity between the meteoroids 
and 2003 EHi. Having established an approximate age from 
backward integrations we then attempt to simulate the forma¬ 
tion of the core of the stream forward in time using the forma¬ 
tion epoch found from backward integration and compare with 
the characteristics of the stream. However, our intention is not 
to provide a complete and detailed picture of all physical char¬ 
acteristics of the stream, rather we aim to demonstrate whether 
the observed overall characteristics of the core of the stream 
can be explained by assuming a relatively recent (a few hun¬ 
dred years) formation age derived from backward integrations 
of individual meteoroids. 


As a test of reliability of our backward integration estimate 
of the age and formation mode of the stream, we compare the 
following theoretical and observed characteristics of the stream; 

1. The timing of the appearance of the stream on the sky 
(around 1835 AD). 

2. The mean position and spread of the geocentric radiant of 
the stream. 

3. The position of the peak of the activity profile of the core. 

4. The width of the activity profile of the core (EWHM 0.6 
days). 


Throughout this work, we use an approach similar to that 
of Gustafson! ( 19891) . That author integrated backward in time 
the orbits of 20 high precision Geminids, along with the parent 
3200 Phaethon, and compared the epochs at which the orbits 
of the Geminids and that of Phaethon intersected. Moreover, 
he calculated the probable meteoroid ejection speed and loca¬ 
tion on the orbit of the parent, and concluded that the Geminids 
are consistent with cometary sublimation that might have taken 
place on Phaethon around 600 -2000 y ears ago. Tor an e xhaus- 
tive description of the method see also Adolfssonl ( 19961) . 
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2. Asteroid 2003 EHi and the Quadrantids meteoroid 
streams 


Asteroid 196256 (2003 EHi) was discovered in March 2003 
by the Lowell Observatory Near-Earth Object Search (LO- 
NEOS) and has been designated as a NEO Amor type aster¬ 
oid. Currently, 2003 EHi moves on a highly inclined cometary- 
like orbit with a Tisserand parameter with respect to Jupiter of 
T 7 = 2.06 but shows no evidence of a cometary activity. 

Recent photometric observations of the asteroid 2003 
EHi give an absolute visual magnitude H - 16.2 mag, 

given by the NASA’s JPL Small Body Database Browser 
(http: //ssd. jpl.nasa. gov/sbdb. cgi). The diameter c an 
be found using the expression given bv ICheslev et al. (2002): 


D{km) = 1329A~‘’'5 x 10 


-0.2H 


( 1 ) 


where A is the geometric albedo and H is the absolute visual 
magnitude of the asteroid, respectively. A typical albedo value 
for asteroids (depe nding on the spectral class of the asteroid) 
is 0.04 < A < 0.4 ( Harris , 1989h . which using Eq. [1] yields a 
diameter of 1.2 km < D < 3.8 km. 

The present orbit of the asteroid 2003 EHi is presented in 
Fig.m along with the orbits of the Earth, Mars and Jupiter 
and the mean Quadrantid orbit. The orbital elements of 2003 
EHi and the mean orbit of the Quadrantids are given in Table [1] 
along with other previously suggested parents, C/1490 Y1 and 
96P/Machholz. The orbital elements of the mean Quadrantid 
stream are taken from .Tenniskens et al.l ( 1997 ). 



X, (AU) 


The observed visual activity profile of the Quadrantid me¬ 
teor shower is presented in Fig. |2] It is a composite (aver¬ 
age) of activity profiles, as deduced from visual observations of 
Quadrantids, in the y ears of 1986, 1987, 1989, 1990 and 1992 
( Rendtel et aU 1993). The core of the activity profile is easily 
seen to be less than a day wide, with a FWHM 0.5 days. 
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Figure 2: The average visual activity profile of the Quadrantid meteor 
shower. The profile is a sta ck of a few favorable observations between 1986 
- 1992 ^Rendtel et al.Lll993h . The peak of the activity is centered around Aq = 
283.2° for equinox of J2000.0. 


Fig. |3] shows the average radar activity profile of the Quad¬ 
rantids (milligram size particles), as observed by the Cana¬ 
dian Meteor Orbit Radar “CMOR” (see Section|4]for details). 
The peak of the activity, located at T© = 283.1° for equinox of 
J2000.0, occurs slightly before the visual meteors peak, but is 
in agreement within uncertainty of the two (averaged) profiles. 
The FWHM corresponds to approximately 0.8 days, slightly 
wider than the visual activity profile (see Fig.|2]). 



K’ deg 


Figur e 1: Orbits of asteroid 2003 EHi (blade line) and the mean Quadrantid 
stream [jenniskens et all 1 19971) (magenta line), viewed from above the ecliptic 
plane. The orbit of Earth is indicated with blue. Mars with red and Jupiter 
with orange colors, respectively. The portions of the orbits, below the ecliptic 
are denoted with a dashed line. The ascending nodes of the orbits of both, the 
mean Quadrantids and 2003 EHi are located close to the Jupiter’s orbit. 


Figure 3: The average radar activity profile of the Quadrantid meteor shower 
as deduced by the Canadian Meteor Orbit Radar “CMOR“. The profile is an 
average profile of all good quality observations between 2002 - 2014. 
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As shown in Fig. [T] the stream of particles producing the 



























Quadrantids intersect the orbit of the Earth at their descending 
node, whereas the ascending node is located near the orbit of 
Jupiter. We emphasize that the mean Quadrantid orbit, as ob¬ 
served from the Earth, does not necessarily represent the orbit 
of the entire physical stream, but rather only that part that in¬ 
tersects the Earth. Meteoroid streams are usually dispersed and 
not all members intersect the Earth, so it is not possible to know 
the orbits of that portion of the stream that does not physically 
intersect the Earth’s path. Therefore we refer to the mean Quad¬ 
rantid stream as the mean of that portion that interacts with the 
Earth’s atmosphere. 

It is also evident that present orbit of the Quadrantid mete¬ 
oroid stream and that of asteroid 2003 EHi are strikingly sim¬ 
ilar. However, the proximity of the parent and stream orbits is 
not a sufficient, nor necessary condition for a relation between 
the two. Eor example, initially close orbits (e.g. stream and 
parent) may follow dramatically different dynamical evolution 
due to the differential gravitational perturbations and solar ra¬ 
diation forces. Conversely, two completely different orbits may 
end up into similar ones due to the aforementioned effect and 
result in a ’’false” sibling relationship between the two orbits. 
Thus, robust models of dynamical evolution of streams must be 
e mployed. 


JenniskensI (120041) were the first to notice the similarity be¬ 


tween 2003 EHi and the Quadrantids and proposed the former 
as the parent of the stream. However, the lack of cometary ac¬ 
tivity on 2003 EHi and the nature of its orbit make it a strong 
candidate for a recently dormant or defunct comet. The real na¬ 
ture, whether cometary or asteroidal, is unknown due to the lack 
of favorable returns and observations, rendering the formation 
mechanism of the Quadrantids uncertain. 

The nature of 2003 EHi will determine the underlying mech¬ 
anism by which it produced the Quadrantids. While cometary 
activity is driven by the sublimation of water ice, asteroi ds may 
also s hed material by other processes. Eor example, Ijewitl 
(2012 h argued that the mass loss mechanism from the surface 
of 3200 Phaethon, another seemingly dormant or extinct comet, 
producing the Geminids meteor shower, is due to solar radiation 
pressure sweeping. 

The inactivity of 2003 EHi points perhaps to an asteroidal 
nature. However, recent spectroscopic and double-station ob¬ 
servations of 51 Quadrantids, revealed that unlike the Geminids 


the Q uadrantids are not as depleted in volatiles (iKoten et al 


20061) . In fact, the authors argued that the Quadrantids fall in 
a category between meteoroids of cometary and asteroid ori¬ 
gin, and concluded that the parent 2003 EHi must be a dormant 
comet. Given 2003 EHi’s comet-like orbit, we will here assume 
that material is released from it by traditional cometary activity 
in our forward modeling. 

3. Meteoroid ejection model 

Besides meteoroid ejection due to cometary sublimation, 
there are a few other potential mechan ism of mass shedding 
from the surface of small bodies (e.g. ( .lewitj. 2012h ') such as 
impact ejection, rotational instability, electrostatic forces, ther¬ 
mal fracture and radiation pressure sweeping. However, in this 


work we limit ourselves to cometary sublimation as the most 
likely mechanisms of formation of the narrow core of the Quad¬ 
rantids. 

3.1. Cometary sublimation 

Presently, the meteoroid production mechanism due to 
cometary volatiles s ublim a tion is relatively well understood, 
see e.g. (Whipple. 19.5 it Jones! 19^ Crifo and Rodiono\l 


1997 


Ma et al. . 2002^ iHughesI bOOOh . All meteoroid ejection 


models to a great extent share the same physical concepts al¬ 
though with slight modifications. 

Perhaps the very first meteoroid ejection model came wit h 
the pioneering work done by Ered Whipple dWhipplel Il9-5lh . 
when he proposed his ”icy-conglomerate” (’’dirty snowball”) 
comet nucleus model. In his model, the comet nucleus consists 
of frozen volatiles - ices with solid refractory particles embed¬ 
ded within the ice. Once the comet n ucleus is close enough to 
the Sun, e.g 3 AU ( Delsemt^ 1982 ). the cometary ices begin 
to sublimate, releasing and dragging along the embedded dust 
particles. In Whipple’s model, the meteoroids leave the surface 
of the comet with speeds: 


Vej =25.4r 




( 2 ) 


where the r is the heliocentric distance in AU, p is the den¬ 
sity of the meteoroid in kg/m^, m is the mass of the meteoroid 
in kg and Rc is the radius of the comet in km. However, o ne ma¬ 


j or dis advantage of Whipple’s model, later improved by iJones 


(Il99-5h . is the assumption of b lackbody lim ited temperature of 
the comet nucleus. Moreover, Jones ( 19951) considered an adi¬ 
abatic expansion of the escaping gas, i.e the sublimation of the 
ices results in cooling of the comet’s surface. The terminal 
speed of the particl es, leaving the surface of the comet, accord¬ 
ing to ljonesi (Il995h is: 


Vej = 29.1 (m/s) 


(3) 


where the variables and the units are the same as in Eql2] 
One major shortcoming of Jones’ model is the assumption 
that the entire nucleus is active. Direct imaging of the nucleus 
of comet 1 P/Halley, by Giotto and Vega comet probes, showed 
that the sublimation of the nucleus is confined to a small frac¬ 
tion of the cometary surface. Nonetheless, a close inspection 
of Eq. |3] shows that the latter expression is very similar to that 
derived by Whipple. This indicates that accounting for the adia¬ 
batic expansion of the escaping gases and eliminating the black- 
body nucleus limitation yields little difference in the final result 
for the terminal velocity of the ejected meteoroids. 

More recently, in their work of mod eling the dynamical evo- 
lution of the Perseid meteoroid stream. lBrown and Jones! (Il998h 
used a slightly modified version of the Jones ( 19951) ejection 
speed as a function of the heliocentric distance. Moreover, the 
probability of particles having ejection speed P(V- Vej) was as¬ 
sumed to be a parabolic distribution, where the meteoroid ejec¬ 
tion directions are distributed isotropically on the sunl it hemi¬ 
sphere . Regarding the meteoroid production, we follow lKresak 
( 19761) . where the production rate is uniform in true anomaly 
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Object 

Epoch (UT) 

a(AU) 

e 

q(AU) 

Q(AU) 

i(deg) 

oj (deg) 

a (deg) 

2003 EHi 

04 Nov. 2013 

3.122 

0.619 

1.189 

5.054 

70.876 

171.354 

282.962 

96P/Machholz 

06 Sep. 2013 

3.034 

0.959 

0.124 

5.944 

58.312 

14.757 

94.323 

1490 Y1 

26 Dec. 1490 

— 

1 

0.737 

— 

51.65 

129.84 

295.89 

Quadrantids 

— 

3.14 

0.688 

0.98 

5.3 

71.5 

171.2 

283.3 


Table ll The osculating orbital elements of asteroid 2003 EHi, comet 96P/Machholz, comet 1490 Y1 and the mean Quadrantid orb it (J2000.0), as g i ven in the 
NASA’s JPL Horizons System (http://ssd.jpl.nasa.gov/sbdb.cgi I. The orbital elements of the mean stream are taken from IJenniskens et al.L|l997h . Q 
and q in the table stand for the aphelion and perihelion distances respectively. 


(v). According to [Brown and JonesI (I1998h , the terminal speed 
at which the meteoroids leave the surface of the comet is given 
by: 


Vej = 10.2 (m/s) 


(4) 


with ejection probability distribution: 


P(V - Vei) = 1 




with 0 < y < 2 Vej and 0 outside 

(5) 

where P(V - Vej) is the probability of finding a meteoroid 
with ejection speed V. 

In the beginning of 2y' century, several interesting works ap¬ 
peared in which authors attempted to infer the ejection speeds 
of meteoroids causing a meteor outburst at the Earth. Observa¬ 
tions of the profile of the shower, given the meteoroids were re¬ 
leased at a certain perihelion passage of the parent, can be used 
to analytically solve for the ejection conditions from the parent. 
The method has ex tensively been applied to the 1999 Leonid 
outburst (see e.g. ( Gockel and Jehn , 20001 Ma and WilliamsL 


2001 : Muller et al. . 200 ih l. This approach, however, has its 


limitations i.e the meteoroids producing the outburst must be 
large (mm - cm) and originat e from relative ly young meteoroid 
trails (a few orbital periods) (lAsherl 120081) . The ejection con¬ 
ditions for much older meteoroids is more difficult (almost im¬ 
possible) to be deduced as planetary perturbations, over several 
orbital periods, become significant thus rendering the recovery 
of ejection conditions impossible. This method, however, is not 
suitable for our simulations for a few reasons: first the most 
probable age of the Quadrantids, based on our backward inte¬ 
grations, (200-300 years 35-55 orbital revolutions of 2003 
EHi) renders it impossible to apply the above method. Sec¬ 
ondly, we do not have a priori knowledge on the physical prop¬ 
erties of the nucleus, as the cometary activity may vary signif¬ 
icantly from comet to comet. This method, however, is worth 
for future investigation. 

Throughout this work, for our forward modeling, we 
choose to use the abov e (Eq 0] and Eq |5]) ejection model by 
Brown and JonesI ( 1998h . We assume that 2003 EHi has a mean 
density of Pr = 800 kg/ m^ - a typical value for a comet nucleus 


(IWeissman et al.L 120041) . while for the radius of the parent body 


we use an average value of 7?^ = 1 km, see Section]^ The me¬ 
teoroids, on the other ha nd, are modeled as sp herical grains of 
.3 


i na, are moaelea as sp nericai grams ot 
density p - 1900 kg/m^ ( Babadzhanovl 120021) . where only the 
radius of the particles is allowed to vary (Section [53]). Due to 


the present lack of cometary activity on 2003 EHi, we assume 
it to be relatively depleted of volatiles such as: CO, CO 2 , CH 4 , 
NHj etc. and thus comprising mostly of water ice (if any) and 
embedded refractory material (meteoroids). The water ice be- 
gins to sublimate a t a heliocentric distance of roughly r = 3 AU 
( Delsemma 1982h . thus the meteoroids, in our simulations, will 
be released from 2003 EHi on an arc of the orbit within 3 AU 
from the Sun, centered at the perihelion of 2003 EHi. Eor ex¬ 
ample, assuming an orbit similar to th at of 2003 EHi, see Tab[p 
and meteoroid ejection model that of lBrown and .Tones! (Il998h . 
the magnitude of the ejection speed of a meteoroid, with bulk 
density of p = 1900 kg.m^^ and mass m - 7.6 x lO^^kg (cor¬ 
responding to a radius of 1 mm), will be Vej ~ 18 m/s at r- 1 
AU, Vej ~ 8.8 m/s at r= 2 AU and Vej ~ 5.8 m/s at r- 3 AU. 


4. Observational Data 


In our study, we use the orbital data of eight high-precision 
Quadrantids, photographed from multiple statio ns by the Czech 
part of the Europea n Fireball Network (“EN”), ( Spurn4 1994 : 
Spurnv et al.L 2007h . In addition to these bright photographic 
Quadrantids, we also use five core radar Quadrant ids observed 
by th e Canadian Meteor Orbit Radar (CMOR) ( .Tones et ^ 
20051) . which sample the smaller members of the stream (mil¬ 
ligram size). These Quadrantids were backward integrated and 
their orbit compared with that of 2003 EHi to search for epochs 
of close encounters which might indicate when the material 
forming the young core of the stream was released. 


4.1. Photographic Quadrantids 

Live out of the eight photographic Quadrantids pertain di¬ 
rectly to the narrow portion of the stream, as they have been ob¬ 
served within a degree from the peak (Aq = 283.2°), which we 
refer to as “core Quadrantids”. The remaining three lie slightly 
outside of the core of the stream with two extreme cases being 
Dec. 31 («i 3° from the peak) and Jan. 8 («; 4.2°) (see Ta- 
ble|2]). Meteors observed a few degrees away from the peak of 
the Quadrantid shower are referred to as “non-core Quadran¬ 
tids”. 

These bright meteors are of an exceptional importance in our 
simulations, as their photometric mass (gram-kilogram) implies 
that the effect of the solar radiation pressure on these meteors, 
given the timescale of our simulations of a few hundred years, 
would be negligible. However, the mass of each individual pho¬ 
tographic meteoroid can be uncertain by a factor of 3 due to our 
poor knowledge of the value of the luminous efficiency t, used 
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for the estimation of the ablated mass. Thus, in order to ensure 
that the effect of the solar radiation pressure has properly been 
taken into, we consider a wide range in the (3 - parameter, where 
P - FrIFc'k the ratio of the solar radiation pressure to the solar 
gravity (Section |5]|. The assumed values for the p - parameter 
are: 0,10“^, 10 10'^. The aim is to encompass the entire pos¬ 

sible range in p, because of the a priori uncertainty in the mass 
of the meteoroids. The magnitude of p can be expressed as: 


P^^^ 3.54 X (“cgs“ units) 

P^^^ 3.54 X 10-" (“Sr units) 

mj 


(6) 


where, for instance, a spherical meteoroid of mass 1 gram 
and bulk density of p = 1.9 gjcrc? has p ^ 2.3 x 10“^. It is ev¬ 
ident that the larger meteoroids are less susceptible to the solar 
radiation pressure and in fact, for bolides (kilogram size) p can 
safely be assumed to be zero. 

Except the photometric mass, the quality of the data for 
each individual photographic Quadrantid is exceptionally good, 
with the uncertainty in the speed (before the atmospheric entry) 
ranging between a; 0.05-0.5 %, whereas the radiant position is 
known to a; 0.02°. Unlike other authors who use the mean or¬ 
bital elements of the Quadrantids as the starting point for in¬ 
tegrations, we use these high precision observable and mea¬ 
surable quantities i.e the geographic coordinates and speed of 
a point on the meteor trajectory and the radiant coordinates. 
These quantities can directly be translated into geocentric state 
vectors for each observed Quadrantid providing a high preci¬ 
sion initial osculating orbital elements used in our backward 
simulations. 


physical information about the meteoroid, a multi-station de¬ 
tection allows for the determination of a meteoroid’s heliocen¬ 


tric or b it (for details, see e . glWebster et al.l (l2004h : I Jones et al 
( 2005h : Wervk and Brown ( 2013 )'). 

During the 2013 campaign, CMOR measured 1440 radar 
Quadrantids orbits (with masses of the order 10“^ - 10 
grams). Out of that list, we have selected 31 Quadrantids with 
heights above 100 km, to minimize atmospheric deceleration, 
and fractional error in the semi-major axis dafa < O.l. Of 
the 31 preliminary selected Quadrantids, we extracted only 5 
with highest signal-to-noise ratio echoes, as well as the least 
uncertainty in the atmospheric velocity 1% where all un- 
c ertainties were found usi ng a Monte Carlo routine approach 
(( Wervk and Brown . 2012t) '). Table[3]lists the observed geocen¬ 
tric quantities and radiants for these high-quality radar orbits. 

While the velocity is known to 1% in our data set, the mass 
estimate of the meteoroid is more uncertain. The ablated mass 
is a function of the speed of the meteoroid V in th e atmosphere 
and the electron line density q ( Vernianil 1973h . created by 
collisions of the evaporated meteoric atoms with atmospheric 
molecules. The major uncertainty in the mass arises from the 
uncertainty in the electron line dens ity q, which translates int o 
an error in the mass of a factor of 3 ( Wervk and Brown . 2012h. 
For a bulk density of a meteor of p =1.9 g/cm^ (IBabadzhanov , 


2002h and mean value for the mass of radar Quadrantids m 


3 X 10“^ grams, see Tabled this yields for a typical p - parame¬ 
ter for a radar meteor p ^ 1.6x10“"^, if we use for the scattering 
efficiency Qpr - 1. Even if the mass for a radar Quadrantid was 
uncertain by a factor of 3, that would translate in error for the p 
- parameter p ^ 2x 10 and p ^ x 10 for an under and over¬ 
estimated mass, respectively. These last results confirm that we 
can safely use values for 10'^< p <10~^, in order to encompass 
the entire possible mass range of the radar Quadrantids. 


4.2. Radar Quadrantids 


5 . Numerical Simulations 


Radar observations of the Quadrantid meteor shower 
have been reported by s everal authors since 1947 
(IHawkins and Almondl Il952h: with some majo r stud ¬ 


ies carried out by e .g. _ IHawkins and Almondl (Il952h : 


Millman and McKinley ( 1953 ): Poole et al. Jl972ir Radar 


observations can be conducted even during daytime and over¬ 
cast weather, and can detect meteoroids of much smaller sizes 
(milligram size) than those responsible for the visual meteors 
(gram-kilogram). Radar observations of the Quadrantids 
complement visual measurements as the shower peaks in early 
January when weather conditions are often poor in the northern 
hemisphere. 

CMOR is a triple-frequency multi-station meteor radar sys¬ 
tem (17.45, 29.85, 38.15 MHz), with the main station located 
near Tavistock, Ontario, Canada (43.264 N, 80.772 W), and five 
other remote sites, recording ~ 4000 orbits per day. Using the 
time delay for common echoes between stations and echo di¬ 
rections of the received signal at the stations, a meteor’s speed 
and atmospheric trajectory can be determined. While single 
station operation registers echoes (where only echo direction 
and ranges from the main site are measured) and provides some 


We explicitly assume that the parent body of the Quadrantids 
is the asteroid 2003 EHi. For each individual Quadrantid, (see 
Table 12 and [3]), we create lO'* hypothetical clones with orbits 
similar to each individual Quadrantid (see Section I5TT) . We 
then numerically integrate the equations of motion of all clones, 
along with the parent, backwards in time (Phase 1). Once the 
most probable formation epoch has been established, via these 
back integrations, then that time window is used for meteoroid 
ejection and forward in time integration of test particle orbits 
until the present (Phase 2). 

5.7. The “clones" 

The key point is to integrate the equations of motion of a 
large number of hypothetical meteoroids (clones) of a Quad¬ 
rantid observed in the present epoch, along with the presumed 
parent body, backwards in time and to statistically attempt to 
determine the epoch of minimum distance between their orbits. 
This approach is more reliable rather than trying to locate the 
nearest epoch, back in time, when the parent and the meteoroid 
actually physically intersect. The latter event is very unlikely to 
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Meteor Number Date A (p h nio R.A Dec. Voo a e i co Q. 

(UT) (°) (°) (km) (g) (°) (°) (km/s) (AU)_(°) (°) (°) 


(J2000.0) 


Core 

Quadrantids 


EN030109 

2009-Jan-3 

11.1406 

48.95368 

93.29 

1200 

227.38 

52.03 

42.57 

2.86 

0.656 

72.18 

176.66 

283.3930 


18:19:14.461 

±0.0003 

±0.00005 

±0.02 

- 

±0.04 

±0.01 

±0.02 

±0.01 

±0.002 

±0.03 

±0.06 

- 

EN040180 

1980-Jan-4 

16.05305 

50.12508 

93.30 

1.9 

232.28 

50.40 

41.92 

2.79 

0.650 

70.98 

168.51 

283.7938 


17:07:34 

±0.00004 

±0.00001 

±0.01 

- 

±0.02 

±0.01 

±0.10 

±0.05 

±0.005 

±0.10 

±0.08 

- 

EN030189 

1989-Jan-3 

14.9824 

48.6881 

93.53 

2.3 

235.90 

49.67 

40.3 

2.82 

0.656 

66.6 

164.6 

283.7938 


4:27:10 

±0.0001 

±0.0001 

±0.01 

- 

±0.07 

±0.03 

±0.3 

±0.15 

±0.02 

±0.3 

±0.3 

- 

EN040192 

1992-Jan-4 

15.98771 

49.92328 

91.33 

80 

227.38 

49.96 

42.54 

2.71 

0.637 

71.87 

174.45 

283.1019 


2:48:30 

±0.00002 

±0.00001 

±0.01 

- 

±0.02 

±0.01 

±0.07 

±0.03 

±0.004 

±0.08 

±0.04 

- 

EN040in 

2011-Jan-4 

15.98771 

49.92328 

91.33 

70 

232.53 

49.67 

41.1 

2.21 

0.642 

68.7 

167.7 

283.2498 


3:19:43 

±0.00002 

±0.00001 

±0.01 

— 

±0.03 

±0.01 

±0.3 

±0.14 

±0.018 

±0.3 

±0.2 

— 

Non-core 

Quadrantids 














EN3n210 

2011-Dec-31 

14.39278 

47.37932 

93.71 

650 

226.77 

51.11 

42.32 

2.88 

0.659 

70.0 

173.10 

279.1733 


3:21:49.71 

±0.00040 

±0.00023 

±0.03 

- 

±0.14 

±0.02 

±0.17 

±0.09 

±0.011 

±0.2 

±0.18 

- 

EN060111 

2011-Jan-6 

16.51110 

48.29288 

80.18 

30 

227.02 

49.86 

42.6 

2.47 

0.602 

71.9 

178.72 

285.3437 


4:36:52.65 

±0.00035 

±0.00025 

±0.05 

- 

±0.08 

±0.08 

±0.3 

±0.12 

±0.019 

±0.3 

±0.18 

- 

EN080179 

1979-Jan-8 

15.61535 

48.96701 

98.02 

60 

233.53 

48.55 

43.21 

3.73 

0.738 

71.61 

170.41 

287.5665 


3:43:55 

±0.00007 

±0.00005 

±0.02 

- 

±0.03 

±0.01 

±0.18 

±0.17 

±0.012 

±0.19 

±0.10 

- 


Table 2l Data for eight high precision Quadrantids, photographed from multiple stations by the Czech part of the European Meteor Network. The table is divided 
into two parts, core Quadrantids - observed within a day of the peak of the shower and non-core Quadrantids - observed within a few days of the peak. The 
first column is the meteor identification number. The second column is the date and time of observation in (UT). The third, fourth and fifth columns are the the 
geographic position of a point on the meteor trajectory (geographic longitude A, geographic latitude <p and altitude h with respect to the sea level). The mass of the 
meteoroid is given in the sixth column. The seventh and eighth columns are the geocentric radiant position for (J2000.0) - Right Ascension a and Declination 6, 
followed up by the pre-atmospheric speed Voo of the meteoroid (ninth column). The last five columns correspond to the heliocentric orbital elements (J2000.0) - 
semi-major axis («), eccentricity (e), inclination (i ), argument of perihelion (co), and longitude of ascending node (Q). All uncertainties are one cr. 


Date 

Time 

(UT) 

A 

n 

<!> 

n 

h 

(km) 

nio 

X 10-3 (g) 

d© 

D 

a 

n 

S 

n 

Va, a 

(km/s) (AU) 

(J2000.0) 

e 

n 

i 

n 

CO 

n 

20130103 

09:11:22 

-82.23 

43.22 

100.863 

1.1 

282.98 

224.2 

46.9 

42.2 

1.81 

0.466 

74.44 

173.4 










±0.7 

±0.07 

±0.023 

±0.41 

±1.7 

20130103 

20:08:57 

-79.47 

43.55 

101.482 

1.7 

283.44 

230.5 

48.4 

40.2 

2.01 

0.512 

69.84 

168.1 










±0.8 

±0.04 

±0.011 

±0.15 

±0.5 

20130103 

20:26:56 

-80.47 

42.83 

100.678 

0.6 

283.45 

230.1 

48.7 

40.9 

2.22 

0.559 

70.56 

169.8 










±0.1 

±0.06 

±0.012 

±0.2 

±0.3 

20130103 

20:50:27 

-78.52 

44.25 

103.118 

5.5 

283.47 

224.1 

46.8 

41.5 

1.87 

0.476 

73.43 

171.6 










±0.1 

±0.11 

±0.031 

±0.59 

±2.1 

20130104 

06:16:11 

-80.95 

42.91 

103.046 

5.7 

283.87 

235.1 

53.6 

39.7 

1.91 

0.488 

68.23 

171.2 










±0.3 

±0.07 

±0.021 

±0.37 

±0.4 


Table 3l Data for five high quality radar Quadrantids, belonging to the narrow core of the stream,observed by CMOR during 2013. The columns represent the 
date of observation, time, geographic longitude A, latitude (p and height h above the ground of point on the meteor trajectory, mass mo of the meteoroid, the solar 
longitude at the time of observation d©, geocentric right ascension a and geocentric declination 6 of the radiant position and the speed Voo of the meteoroid above 
the Earth’s atmosphere. The last four columns are the computed heliocentric orbital elements of the meteoroid, semi-major axis a, eccentricity e, inclination i and 
argument of perihelion co. The average uncertainty in the radiant position is ± 0.2°. 


7 







occur, due to the inherent uncertainties in the position and ve¬ 
locity of the bodies, as well as computational errors (truncation 
and round-off) throughout numerical integrations. Even though 
the two bodies are seen to intersect in the simulations, that does 
not necessarily imply that we have identified the epoch of mete¬ 
oroid release, due to the reasons mentioned above. A simple ex¬ 
ample can be used to demonstrate that: Even one has calculated 
the semi-major axis a of a meteoroid orbit with an uncertainty 
of, e.g a ± da - 2.855 + 0.013 AU (as in the case of meteor 
EN030109 Table |2]), over several revolutions of the meteoroid 
about the Sun, the uncertainty in its position will dramatically 
increase. Using and differentiating Kepler’s third law, it can be 
demonstrated that; 


dT = 


3a^^^da 

2vr^’ 


(7) 


where dT is the absolute uncertainty of the period of an 
orbit around the Sun. This yields a relative uncertainty of 
dT jT ss 7 X 10“^ even for 13 — 0, corresponding to dT » 12 
days per period or one orbit within » 700 years, comparable to 
the time scales considered here. It is noteworthy that absolute 
uncertainty in dT increases for finite values of p. Thus, the data 
is not precise enough to determine when the meteoroid and the 
parent actually intersect: such studies will have to await even 
higher quality measurements. Nonetheless, we can use the in¬ 
tersection of the meteoroid and parent orbits to examine when 
such ejection could have occurred, as this is a necessary condi¬ 
tion for meteoroid stream production. 

Eor our study, for each observed Quadrantid we create 10^ 
hypothetical clones, randomly selected from a six-dimensional 
Gaussian distribution of the six directly observed and mea¬ 
sured quantities: geocentric radiant coordinates (a, 6), the ge¬ 
ographic longitude, latitude and altitude (A, <p and h) and the 
observed velocity at a point on the meteor trajectory (Vots)- 
These values can be translated directly into the initial helio¬ 
centric state vectors {Vu, Viy, Vi^] and {Xi, T/,Z,} of each clone 
and are used as initial conditions for our backwards numeri¬ 
cal integrations of the equations of motion. Each dimension of 
the six-dimensional Gaussian, is centered around the nominal 
values with one standard deviation 1 cr being equal to the mag¬ 
nitudes of the uncertainties of the initial variables (see Table|2l. 
However, at the time of our simulations, we do not have infor¬ 
mation about the covariance between the errors in our initial 
condition quantities, thus we assume that the errors are uncor¬ 
related, which may not be the case. 


5.2. Phase 1: Backward integrations 

The equations of motion of each clone along with the par¬ 
ent are then integrated backwards in time for 1000 years. The 
length of the backward integrations of 1000 years is chosen on 
the basis of the assumption of relatively young age of the central 


portion of the Quadrantid meteor oid stream (iJenniskens et al. 


1997 : Wiegert and Brown . 2005h and because we found that 


similar Quadrantid orbits in the current epoch integrated back 
in time, do not start to significantly disperse until s; 1300 AD. 
In addition, we aimed to check whether the age of the core of 


the Quadrantids could be as old as 5 00 - 1000 years as pre 


viously argued by some authors e.g. (IJenniskens et al.L 11997 


Williams and Wu . 1993h . We note that integration for 1000 


years back in time is beyond the Lyapunov time, however the 
mean orbital elements of the clones of 2003 EHi, at the epochs 
of 1000 AD and 1500 AD, are used for forward integrations in 
order to compare the theoretical characteristics of core of the 
Quadrantid stream against the observed ones. Effectively, we 
explore whether the observed characteristics of the core of the 
Quadrantids can be explained by meteoroid ejections from 2003 
EHi circa 1000 AD or 1500 AD. However, the results from for¬ 
ward integrations should be treated with caution, in particular 
when the initial orbital elements of the parent have been se¬ 
lected from backward integration beyond the Lyapunov time. 
Therefore, the results must be treated from a statistical point of 
view. 

Throughout the integrations, we consider the perturbations 
on the meteoroids and 2003 EHi from each planet and also ac¬ 
count for the planets’ mutual interaction. We used the JPL’s 
”DE405” version of the planetary ephemeris for generating the 
initial positions of planets. The mutual interactions among the 
meteoroids and 2003 EHi are neglected, i.e they are considered 
as test particles. Thr oughout the inte grations we use Everhart’s 
RADAU algorithm ( Everhart . 1985h . During the first 30 days 
of the integration of the orbits of the clones, i.e. when the mete¬ 
oroid is in the Earth’s vicinity (several Hill radii, where the Hill 
radius of the Earth ss 0.01 AU), we use a fixed time step of 1 
minute, with the gravitational influence of the Moon being sep¬ 
arately taken into account. When the meteoroid is sufficiently 
far away from the Earth, we then increase the time step from 1 
minute to 1 day (from -30 days to -1000 years) in order to speed 
up the integrations. The orbital elements for each hypothetical 
clone and 2003 EHi are output at every 10 years. Eurthermore, 
for all meteoroids, except for EN30109 (the largest at » 1.2 
kg.), we examine four different values for the ratio of the solar 
radiation pressure to the solar gravity ”/J - parameter”, P - 0, 
P = 10-3,/? = 10-4 and/? = 10-^. 

Throughout the backward integrations, we compute the mini¬ 
mum distance between the orbits (MOID) of each clone and the 
parent 2003 EHi. The aim is to statistically determine the epoch 
when the MOID between each observed Quadrantid and aster¬ 
oid 2003 EHi was at a minimum. We assume that the spread- 
ing time in the m ean anomaly («; 100 years for JFCs and NEOs 
( Tancredil 1998 )) of meteoroids on their orbit is much shorter 
than the spreading time in the other orbital elements. Back¬ 
wards integrations within a few Lyapunov times can reliably 
provide information on the minimum distance between the or¬ 
bit of the parent and each Quadrantid, though longer ones may 
not. Our calculated Lyapunov time for asteroid 2003 EHi is s; 
80 years which is consist ent with the values for JFC and NEOs 
found by Tancredil ( 1998 )). 

In addition to the computation of the MOID between the par¬ 
ent and each Quadrant id, we also used the standard orb ital sim- 
ilarity functions Ds h (ISouthworth and Hawkins! Il963h and D' 
( Drummondl 1981 ). between the orbits of the parent and each 
clone as a check on our MOID results and found similar be¬ 
havior. In an attempt to identify a probable formation mode of 
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the core of the Quadrantid stream, we also compute the true 
anomaly 6 at the MOID of the orbit of 2003 EHi and each 
clone. Meteoroid ejection due to cometary ice sublimation must 
be confined within an arc of the orbit when the parent is suffi¬ 
ciently close to the Sun, i.e within r = 3 AU. The backwards in¬ 
tegrations of the equations of motion for the parent, 2003 EHi, 
indicate that its orbital elements do not vary dramatically over 
200 years from their present values. If the narrow portion of 
the Quadrantids was created through cometary meteoroid ejec¬ 
tion around 1800 AD, that must had happened within an arc of 
6 ss + 130° from perihelion which is roughly equivalent to a 
heliocentric distance of r « 3 AU. 

In order to further constrain the meteoroid release mode, we 
also compute the relative velocity between the parent and each 
clone at the MOID. The aim is to compare whether the ejection 
velocities within the arc of water ice sublimation are consistent 
with cometary dust ejection speeds (see Section^. Although, 
all meteoroid ejection models yield slightly different ejection 
velocities (a few tens to few hundred m/s), nonetheless they all 
agree that the latter are unlikely to exceed ~ 1 km/s. Thus, 
if the relative velocity between the parent, 2003 EHi, and each 
observed Quadrantid is within few-hundreds m/s, we may argue 
for cometary origin of the core of the Quadrantid stream. 


5.3. Phase 2: Forward integrations 

Using the results of the backward integrations to provide a 
time window for the likely formation epoch of the core of the 
Quadrantid meteoroid stream, we next use ejection models (de¬ 
scribed below) to generate a hypothetical meteoroid stream and 
integrate the equations of motion of the meteoroids forward in 
time. 

Based on the results of Phase 1, we assume that meteoroid 
ejections took place between 1700 AD and 1900 AD, and we 
use this time window to eject meteoroids from the parent 2003 
EHi between 1700 - 1900 AD. We explicitly assume that par¬ 
ent, 2003 EHi, was an active comet prior to 1900 AD, supply¬ 
ing meteoroids to the Quadrantid meteoroid stream, and that 
activity ceased after the beginning of 20'^ century. The latter 
assumption will be motivated by our backwards integrations as 
described later. Thus, we effectively test the hypothesis whether 
meteoroids ejected during between 1700 - 1900 AD could re¬ 
produce the presently observed narrow structure of the Quad¬ 
rantid stream (see Section|6]for details). 

We hrst test what is the most likely epoch, at which mete¬ 
oroids have to be released from 2003 EHi, so the resulting me¬ 
teoroids start intersecting the Earth’s orbit around 1835 AD i.e 
when the stream was first noticed. Eor this purpose, we eject 
10“* meteoroids, of both radar (lOOjum - 1 mm) and visual sizes 
(1 mm - 1 cm), at various epochs as a single outburst from 
the parent and propagated their equations of motion forward 
in time. The ejection model that we used for the single out burst 
meteoroid ejection is the one by Brown and JonesI ( 1998 ). de¬ 
scribed in Section [3] We note here that the epochs, at which 
radar and visual size meteoroids are released from the parent 
so they can reach the Earth around 1835 AD, are slightly dif¬ 
ferent. This is due to the different orbital evolution of micron 
and centimeter sized meteoroids, for which the the magnitude 


of the solar radiation pressure force is different. Once the ’’cor¬ 
rect” time window of meteoroid release has been identihed (i.e 
the resulting meteoroids reach the Earth around 1835 AD), we 
then use that time window as a starting point for continues me¬ 
teoroid ejection over multiple perihelion returns of 2003 EHi. 
In the case of visual size particles this time window is between 
1780 - 1786 AD, whereas for radar size meteoroids it is between 
1790 - 1796 AD (see Section|g. 

In the case of meteoroid ejections over multiple perihelion 
returns of the 2003 EHi, 10"^ meteoroids are released at each 
perihelion passage of the parent (within an arc of 3 AU, cen¬ 
tered at the perihelion of 2003 EHi), both at radar and at visual 
sizes, to compare with both observation techniques. The me¬ 
teoroid ejection speeds considered are described in Section [3] 
The number of particles, in a given size bin, were chosen uni¬ 
formly in the logarithm of the size. During the integrations, the 
orbital elements, the resultant radiant position and geocentric 
velocity for the meteoroids are computed for each day. That 
allows us to closely examine the activity profile of the resultant 
stream as a function of the meteoroid size, i.e radar and visual 
meteors, and compare it with the observed activity prohle, ra¬ 
diant position and dispersion and geocentric speed. 

We would like to note, that we also performed sample simu¬ 
lations assuming the extreme values for the radius of 2003 EHi, 
i.e Rc -2 km (see Section|2]l as well as the density of the mete¬ 
oroids within 0.8 g/cm^ < p <3 g/cm^. However, the variation 
of the radius of the parent, did not show a noticeable difference 
in the hnal results. Furthermore, we also repeated the simu¬ 
lations with all variables as above, i.e parent’s radius and den¬ 
sity of meteoroids, number of meteoroids ejected per perihelion 
return of 2003 EHi, meteoroid masses and /?-values etc., thus 
replacing only the ejection speed model with the one result¬ 
ing from h ydrodynam i cal study of the cometary circum-n ucleus 
coma by ( Crifo et al. . 1995^ Crifo and Rodionov . 1997 ). The 
latter model is generally known to yield sl ightly lower me ¬ 
teoroid ejection speeds, as compared to e.g IWhipnlel (Il95ll) : 

Effectively, this allowed us to test the 


Brown and Jone: 


s peeds, 

idl998h. 


significance of the magnitude of the ejection speeds on our final 
results, keeping all other parameters same as above. However, 
the overall results did not show a modest difference, thus imply¬ 
ing that the considered ejection model has negligible effect on 
the final results, over the time scale of integrations we are con¬ 
cerned with. We thus, did not see an y strong reason to further 
investigate which ejection model e.g. ( Brown and JonesL 1^98) 


{mi; 

Crifo and RodionovI (Il997h . Therefore, in this work we only 


should be considered ove r another e.g. ICrifo et al. 


present the results from the simulation, c arried out assum¬ 
ing m eteoroid ejection speeds modeled by ( Brown and JonesL 


m^- 


6. Results 

6.1. Phase 1: Meteoroid release epoch 

In Phase 1, we integrated 10^ hypothetical clones back¬ 
wards for each observed Quadrantid (eight photographic and 
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five radar). For brevity we show the results, from backward in¬ 
tegrations, for only three core and two non-core Quadrantids. 
Moreover, we show the results for only y6 = 0, as the other in¬ 
stances of p yielded similar outcomes. The simulations are thus 
almost insensitive to the meteoroids’ mass range, given the time 
scale of our simulations. Despite solar radiation pressure force 
being non-negligible for micrometer and millimeter sized parti¬ 
cles, we do not hnd that modestly larger beta values noticeably 
change the hnal results. 

Fig. |4] shows the evolution of the MOID (left panels) and the 
relative velocity at the MOID (right panels), between the orbits 
of 2003 EHi and 10"^ clones, for three core Quadrantids. The 
y - dimension of each pixel in Fig.|4]is equal to ~ 0.005 AU in 
the left panels and ~ 0.01 km/s in the right panels, whereas the 
x-dimension is 10 years. It is clearly seen, from Fig.lH that the 
MOIDs between the clones of each integrated core-Quadrantid 
and 2003 EHi show a deep minimum 200 - 300 years before 
the present (a; 1700 - 1800 AD), with little dispersion in the 
MOID at that particular epoch. The same general tendency of 
the MOID, with a minimum between 1700 AD -1800 AD, was 
observed even for particles with extreme values of p - 10“^. 

Similarly, Fig. |6] shows the evolution of the MOID (left 
panel) and relative velocity at the MOID (right panel), as above, 
for ”non-core” (observed outside the narrow peak of the stream) 
meteoroids. Unlike core-Quadrantids, their non-core counter¬ 
parts did not show an obvious minimum of the MOID. The 
exception was non-core meteor number EN060111 which did 
show a weak minimum in the MOID around 1650 AD but with a 
relatively high dispersion. We suggest that the non-core Quad¬ 
rantids are much older, though the Lyapunov times are such 
that longer backwards integrations to confirm this hypothesis 
are problematic. 

As an additional test as to the age of the core of the Quadran- 
tid str eam, we also computed th e stand ard similar ity functions 
Dsh (ISouthworth and Hawkinsj Il963h and D’ (iDrummond , 


19811) . between the orbits of the clones of each individually ob¬ 


served core-Quadrantids and 2003 EFIi Fig. |5] However, here 
we present the results for only 3 bolides, EN030109, EN040180 
and EN040192. The results for the rest of the core-Quadrantids 
yielded similar results. The size of each pixel in the left panels 
is (5x = 10 years, 6y - 0.01, whereas in the right panels {D’ 
criterion) (5x = 10 years, 6y - 0.005. Both similarity criteria 
yielded an unambiguous deep minimum between 1700 - 1800 
AD, similar to the MOID criterion, for the core-Quadrantids. 
However, in the case of non-core Quadrantids there was not a 
clear minimum of the similarity function, the results which we 
have omitted here. 

The effect of varying the p - parameter had little effect, that is 
the evolution of both, Dsh and D' criteria, yielded a minimum 
between 1700 - 1800 AD. The latter result is not surprising, 
given the time scale of our simulations and in particular the 
time of interest (200 - 300 years), i.e a few Lyapunov times. 

We present the results for only one radar core-Quadrantid, 
since all other cases yielded similar results. From Fig. |7] it 
can be seen that, as in the case of the photographic core- 
Quadrantids, the minimum value of the MOID between the or¬ 
bits of the parent and the hctitious clones of meteor 20130103- 


20:26:56 is reached around 1700 - 1800 AD. The simulations 
for different p - values also yield similar results. Fig. [8] shows 
the evolution of the orbital similarity functions, Dsh and D', 
between 10^ clones of the radar core Quadrantid 20130103- 
20:26:56 and 2003 EHi. In this case, the minimum of the sim¬ 
ilarity criteria was somewhat less obvious but still around 1800 
AD. We note that radar orbital measurements are of lower pre¬ 
cision compared to photographic or video observations. Never¬ 
theless, the minima of the similarity functions, occurring after 
1700 AD supports the hypothesized young age of the core of 
the Quadrantids. 

As was the case with the photographic core Quadrantids, 
variation of the ratio of the solar radiation pressure to the solar 
gravity p did not change the overall results for radar meteoroids 
as to the position of the minima of both similarity functions. 
The latter result seems to be reasonable, given the time window 
we are interested in (200-300 years). 

All lines of evidence from the backward integrations in Phase 
1 point strongly to an origin in the past 200-300 years for the 
core Quadrantids. 

6.2. Phase 1: Formation mechanism of the Quadrantid mete¬ 
oroid stream 


Another goal of this work is to constrain the formation mode 
of the core of the Quadrantid meteoroid stream. More precisely, 
we aim to test if it is consistent with the hypothesis of cometary 
origin. 

We recall that existing dust ejection models from comet nu¬ 
clei suggest that the ejection velocities are much smaller (few 
tens-hundreds m/s) than the orbital speed of the comet (few tens 


ceeding 1 km/s (see e.g., ( 

Whinole 

1950l 

1951; .lones. 1995 

Crifo et al.l 1995; Crifo and Rodionov, 

1997h). 


In order to test if the Quadrantid meteoroid stream may have 
resulted from a c ometary activity o n 20 03 EH| we follow an ap¬ 
proach similar to Gustafson ( 1989 ) and lAdolfsson ( 1996 ), com¬ 
paring the relative velocity at the MOID between 10"^ clones of 
each photographic and radar core-Quadrantid and 2003 EHi. It 
has to be noted, though, that this approach has its limitations 
and should be used carefully. The uncertainties in all measure¬ 
ments e.g. position of the meteoroids and velocities are hnite 
which over time will tend to increase. Therefore we empha¬ 
size that, we utilize this method assuming the age of the ejected 
meteoroids is less that a few Lyapunov spreading times of the 
meteoroids along their orbits. Secondly, there are uncertainties 
in the calculations of the meteoroids’ initial (pre-atmospheric) 
masses, where the latter come into the equations of motion pri¬ 
marily by their/?- values. Thus, the results from the simulations 
must be treated from a statistical point of view. Nevertheless, 
assuming that the core of the Quadrantid meteoroid stream is 
relatively young , compared to the Lyapunov spreading time, 
we should be able to approximately constrain its age, treating 
the results from the backward and forward integrations in a sta¬ 
tistical sense. 

The results for the relative velocity between 2003 EHi and 
10^ clones of each photographic Quadrantid is presented in 
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Figure 4: The Minimum Intersection Distance (MOID) (left panels) and the relative velocity at the MOID (right panels) between the orbits of asteroid 2003 EHi 
and 10^ clones for each the three photographic meteors, belonging directly to the core of the Quadrantid meteoroid stream. The pixels in the figure are color coded 
in the hue of the blue color, with a darker blue corresponding to a greater number of clones. The red curve in the right panels correspond to the median value of 
the MOID at the given epoch. The relative velocity at the MOID is presented in a logarithmic scale in the panels on the right hand side of the figures. Clones with 
MOID > 0.01 AU are not plotted in the right hand panels. 
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Figure 5l The evolution of the similarity criteria, Dsn ^Southworth and Hawkind.ll963l) (left panel) and D' jPrummonA 1 198 iT) (right panel), for 10"^ clones of 
each Quadrantid as a function of time. The red curve corresponds to the median value of the similarity functions, Dsh and D', respectively. 
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Figure The Minimum Intersection Distance (MOID) (left panels) and the relative velocity at the MOID (right panels) between the orbits of asteroid 2003 EHi 
and 10“^ clones for two photographic meteors, not belonging directly to the narrow core of the Quadrantids. Clones with MOID > 0.01 AU are not plotted in the 
right hand panels. 
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Figure 1\ Evolution of the Minimum Orbit Intersection Distance (MOID), between the orbit of 2003 EHi and ten-thousand fictitious clones of the radar meteor 
- 20130103-20:26:56 over a period of one-thousand years. This radar meteor belongs to the core of the Quadrantids. The red curve represents the median value of 
the MOID at a given epoch. Clones with MOIDs > 0.01 AU are not plotted in the right hand panels. 
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Figure 8l Evolution of the orbital similarity criteria, Dsh and D', between the orbit of 2003 EHi and ten-thousand fictitious clones of the radar meteor - 20130103- 
20:26:56 over a period of one-thousand years. This radai' meteor belongs to the core of the Quadrantids. The red curve represents the median value of the similarity 
functions. 
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Fig.|4](right panels) and radar Quadrantids in Fig.|7](right pan¬ 
els). The epoch of interest is the one centered between -200 and 
- 300 years from the present. Relative velocities are computed 
only for clones with MOlDs less than 0.01 AU from the orbit 
of2003EHi. 

During the minimum of the MOID (left panel), the relative 
velocity at the MOID between the clones of each core Quadran- 
tid and the parent 2003 EHi, range from about 0.1 to 2 km/s. 
These values are on the high side for cometary ejection pro¬ 
cesses, as the ejection speeds of meteoroids leaving the surface 
of 2003 EHi, given the perihelion distance, even for microm¬ 
eter meteoroids are unlikely to exceed 200 m/s. However, the 
process of generating clones within the uncertainties necessar¬ 
ily generates many incorrect orbits; only one clone within the 
ensemble is the real particle. Thus, the existence of even one 
clone with a relative speed at the MOID suggests a possible 
cometary origin. 

The lowest relative velocities within the clones of the core 
Quadrantids and 2003 EHi, range from 200 m/s - 800 m/s, with 
the majority of the clones having relative velocities exceeding 
1 km/s. Only 2 clones of bolide EN040180 attained a very low 
relative speed of a; 80 m/s while most had speeds as great as a 
few km/s. We conclude that these results are marginally con¬ 
sistent with cometary ejection processes but can not be used as 
definitive proof for a cometary sublimation origin for the core 
of the stream. 

Similar results are found for the radar core Quadrantids 
(Fig. 0. Here we see even larger dispersion in the ejection 
speed from the parent, ranging from 0.8 km/s to ~ 5 km/s, 
around 1800 AD. However, only between 1-4 particles demon¬ 
strated relative speeds of 0.8 km/s, while the majority had 
speeds above 1-3 km/s. This result is poorly consistent with 
the expected meteoroid ejection speed from comets, i.e it is un¬ 
likely that meteoroids would leave the surface of the parent with 
speeds of 800 m/s, at a heliocentric distance of 1 AU under 
cometary sublimation alone. We recall, however, that the errors 
in the initial values of the position and speed of the radar mete¬ 
oroids in the Earth’s atmosphere are larger than those obtained 
by photographic observations. 

As an additional constraint, as to the more likely formation 
mechanism of the core of the Quadrantids, we also calculated 
the true anomaly at the MOID, of the orbit of 2003 EHi, for 
10"* clones of each Quadrantid (see Section 15721) . Sublimation 
of cometary ices must take place closer to the Sun (within ~ 3 
AU for 2003 EHi, assuming sublimation of water ice only). 

Eig. l9] and ITOl shows the true anomaly at the MOID for ten- 
thousand clones as a function of time of the photographic me¬ 
teoroid EN040180 (1.9 g) and radar event 20130103-20:26:56, 
both of which are core-Quadrantids. It can be seen that, for 
the epoch 1800 AD, there is a greater probability that the true 
anomaly at the MOID, on the orbit of 2003 EHi, was between 
V a; +60° from perihelion (translating to a heliocentric distance 
r ail.3 AU), which is consistent with cometary activity, result¬ 
ing in meteoroids ejection from 2003 EHi. 

Examination of the results for the rest of the photographic 
and radar meteors in our sample were similar to those in Eig. |9] 
and [To] with some showing lower or higher dispersion in true 



Time (yr), 0 = present 



0 2497 4994 7491 9988 

pixel value 


Eigure 9; The true anomaly of 2003 EHi at the MOID of 10'* fictitious me¬ 
teoroids (clones) of the core-Quadrantid photographic meteor - EN040I80 for 
a time span of 1000 years. 
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Eigure 10; The true anomaly of 2003 EHi at the MOID of 10'* fictitious 
meteoroids (clones) of the radar meteor - 20130103-20:26:56, for a time span 
of 1000 years. 


anomaly at the MOID at the epoch of interest (1700 -1900 AD). 
The effect of varying the /? - value had little effect on the final 
outcome. The highest median value of the true anomaly at the 
MOID was V a! +120° from perihelion (roughly corresponding 
to a heliocentric distance r ^ 2.4 AU). Nonetheless, all results 
showed a higher probability that the true anomaly at the MOID 
on 2003 EHi occurred close to perihelion, consistent with me¬ 
teoroid ejection due to sublimation of ices. 

The location of the MOIDs are certainly consistent with 
cometary ejection processes, while the relative velocities are 
low enough in some cases, but not all. Since, it is not possi¬ 
ble to know the ’’true” backward evolution of any of the ob¬ 
served Quadrantid meteoroids, due to the inherent uncertainties 
in their position and speed, thus if even one clone has relative 
speed consistent with cometary ejection, it may point to origin 
consistent with cometary sublimation. We conclude here that 
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cometary ejection processes may be a likely source of the core 
Quadrantids, or at least that we have sufficient grounds to pro¬ 
ceed with forward modeling of the stream under this assump¬ 
tion. 


6.3. Phase 2: Meteoroid ejection 

In order to test the results obtained by backwards integra¬ 
tions, we employ the ’’direct approach” by ejection of a large 
number of hypothetical meteoroids from the parent, and inte¬ 
grate them forward in time. By the ’’direct approach” we test 
if the observed average physical characteristics of the Quadran- 
tid stream, described in Section [T]and|2] are consistent with 
cometary ejection from 2003 EHi around 1700 - 1900 AD. 

We briefly recall the constraints that we are attempting to 
match by the forward integrations. 

1. The timing o f the appearanc e of the stream on the sky 
(around 1835 AD) dOuetelei . 1839ll . 

2. The mean position and spread of th e geocentric radiant of 
the stream (e.g. djenniskens et al.Lll997h ). 

3. The position of the p eak of the activity profile of the core 
(e.g. ( Rendtel et al. . 1993 0. 

4. The wi dth of the activity profile of the core (FWHM 

0.6 days) (Hughes and Tavlot . 1977t Rendtel et al. . 1993 : 


Jenniskens. 2006). 


We first show the results for 10"^ meteoroids, ejected from 
2003 EHi, at different epochs in a single outburst (at a single 
point on the orbit). The aim is to obtain a time windows as to 
when meteoroids must be ejected so they can reach the Earth’s 
orbit around 1835 (see Section 1531) . Thus, meteoroids ejected 
prior to that time window will reach the Earth before 1835, i.e 
when the shower was first observed. Eurthermore, we have to 
identify the epoch of ejection when meteoroids cease reaching 
the Earth. We then use these time windows for meteoroid ejec¬ 
tion, at every perihelion passage in order to closely inspect the 
characteristics of these synthetic meteor showers. 

Eig. nniupper panels) shows the evolution of the descending 
nodes of visual size meteoroids ejected from 2003 EHi, in a sin¬ 
gle perihelion passage, in 1780 and 1786, respectively. Also, 
Eig. [TT] (lower panels) show the mean radiant position, of the 
simulated meteoroids, in sun-centered ecliptic longitude A- Aq 
and latitude (3. In the figures, only meteoroids having their de¬ 
scending node within 0.01 AU, from the orbit of the Earth are 
plotted so we can sample only those particle which may inter¬ 
sect the Earth. It is evident that meteoroids ejected prior to 1780 
arrive near the Earth’s orbit somewhat before 1835 (see Eig.fTTl- 
a), which is inconsistent with the timing of the shower’s appear¬ 
ance (circa 1835 AD). However, it is difficult to argue when 
exactly meteoroids must had been released from the parent in 
order to arrive at the Earth at the right time (circa 1835). This 
is partly due to the unrealistically low number of the simulated 
particles considered in our simulations, as well as the lack of 
exact conditions during meteoroid ejection. In fact, in Eig.fTTl- 
a, although meteoroids ejected in 1780 arrive at the Earth a few 
decades before the shower’s first appearance in the sky (1835), 
the activity level could had been much lower and the shower 
could had easily gone undetected. Therefore, we choose to use 
as initial ejection epoch the year 1780 instead of 1786. 


The theoretical radiant position of the meteoroids (Eig. [IB is 
plotted along with the mean observed positions and dispersions 
of the Quadrantids as obtained by photographic and video me¬ 
teor surveys. The SOMN Quadrantids were selected as being 
potential Quadrantids if their radiant was within 5 degrees of 
the nominal Quadrantid radiant and if their speed was within 
10% of nominal Quadrantid speed. Similarly, Eig.fTTT-b shows 
that meteoroids ejected around 1786 appear at the Earth just 
slightly after 1835 AD. The slight timing inconsistency is per¬ 
haps due to small number statistics. In fact, out of 10"^ particles 
ejected, only about 30 reach the Earth. Nevertheless, within the 
statistics considered in our simulations, the time of ejection and 
first appearance of the stream on the sky is in a good agreement 
with the minimum of the MOID in the backward integrations 
(Section [5.21) . The theoretical radiant position seems to be in 
a reasonable agreement with the observations by SOMN, al¬ 
though there is slight discrepancy with the DMS photographic 
and video data (Eig.[TTI-a. b lower panel). Thus particles ejected 
from the parent in 1786 by cometary processes produce a rea¬ 
sonable match with the observed shower properties. However, 
the descending nodes of meteoroids ejected after 1882 AD start 
to gradually move outward from the Earth’s orbit, and eventu¬ 
ally stop intersecting the Earth. The reason for that is the sec¬ 
ular precession of the orbital elements, causing the heliocentric 
distance of the descending nodes of the meteoroids to recede 
slowly from the orbit of the Earth. Thus the current visual ac¬ 
tivity of the core of the Quadrantids shower can be explained 
by cometary activity of 2003 EHi between the years 1786 and 
1882. The parent cannot have been active earlier because the 
Quadrantids shower would have become active earlier than is 
recorded in historical records. Whether the parent may have 
been active after 1882 remains unconstrained, except by the 
fact that dust production must have eventually diminished to 
the currently observed inactive state. 


Source 

A - Aq 
(deg) 

yS 

(deg) 

D- (2 - 2o) 

(deg) 

0-/3 

(deg) 

Photographic 

—DMS 

277.8 

63.5 

2.4 

1.1 

Video 

—DMS 

278.7 

63.6 

3.5 

0.9 

—SOMN 

275.8 

63.6 

6.2 

1.7 

Radar 
—CMOR 

277.2 

63.2 

6.8 

3.1 


Table 4; Radiant position and dispersion of tFie Quadrantid meteor shower 
obtained by photographic, video and radar techniques. The coordinates of the 
average radiant are given in the sun-centered ecliptic frame - ecliptic longitude 
(/I — /to) 2 nd ecliptic latitude fS. 


Similarly, Eig. [12] shows the evolution of the descending 
nodes of radar size (100 jim - 1 mm) meteoroids, ejected from 
2003 EHi in 1790 AD and 1796 AD. Meteoroids, ejected prior 
to 1796, appear at the Earth as early as 1800 AD. However, we 
do not really know when radar size particles first reached the 
Earth, as the earliest Quadrantid meteor radar observations date 
back only about six decades. The radiant position Eig. \T7[ a 
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Figure 11 ; Evolution of the descending nodes of meteoroids (upper panels) ejected from 2003 EHi in 1780 and 1786. The red line marks the year of 1835 AD, i.e 






















(lower panel) seems to fit well the observed mean Quadran- 
tid radiant as measured by the CMOR. Thus 1796 AD seems 
to be in a good agreement with the minimum of the MOID, 
between the orbits of radar meteoroids and that of 2003 EHi 
(Section 15.21) . Meteoroids ejected after 1886 AD, from 2003 
EHi, do not presently intersect the Earth’s orbit due to, as in 
the case of visual meteoroids, secular precession of the orbital 
elements of the meteoroids. 

We next present the results of meteoroids ejected over mul¬ 
tiple perihelion returns of 2003 EHi for both, visual and radar 
size particles. According to our analysis from meteoroid ejec¬ 
tions at a single point on the orbit of the parent, we obtained that 
visual size meteoroids appear and intersect the Earth if ejected 
roughly between 1780 - 1882 AD, whereas for radar size mete¬ 
oroids this time window is 1796 - 1886. In this context, approx¬ 
imately 10"^ meteoroids were ejected per perihelion passage of 
the parent, within an arc of 3 AU from the Sun. In the case of 
the visual size particles, the meteoroids were ejected over ~ 18 
perihelion passages of 2003 EHi with total number of ejected 
particles a; 1.8x10^, whereas in the case of radar size mete¬ 
oroids the the equivalent orbital revolutions of 2003 EHi are 
16, totaling in 1.6x10^ ejected meteoroids. 

The stacked theoretical activity profile (blue boxes) of the 
visual Quadrantids ejected between 1780 and 1882 AD, is pre¬ 
sented in Eig[T3] It can be seen that the location of the peak ac¬ 
tivity of our simulated Quadrantid stream, as well as the EWHM 
of the core activity, match fairly well with the observations. 
The wider portion of the observed activity profile though is not 
well reproduced; we suspect these wings are likely much older 
than the timescales (of order a few hundred years) that we are 
concerned with it in this work. The theoretical radiant of me¬ 
teoroids, ejected between 1780 and 1882 AD, shows a fairly 
good match with the observed mean position and dispersion of 
individual visual Quadrantids (Eig. [Hi. 

Fig- US] shows the stacked theoretical and observed radar ac¬ 
tivity profiles. The theoretical width of the core is narrower than 
the observed one, suggesting that there may be an older enve¬ 
lope of meteoroids encompassing the young core. That is, while 
the central portion appears to be relatively young, the wings are 
not. 

Fig. [m shows the simulated distribution of radar size Quad¬ 
rantids ejected between 1796 and 1886. In this case, the theo¬ 
retical radiants yield an even better match to the observed one 
than the visual Quadrantids (Eig. [a. However, meteors de¬ 
tected by radar techniques usually have a greater radiants scatter 
than photographic and video techniques. Thus, the better radi¬ 
ant match may simply be a result of a larger observed radiant 
dispersion. 

Some authors have argued for a much older age (older 
than 1000 years) of the Quadrantid meteoroid stream 
( Williams et al. . 19791: Hughes et al. . 19791 : Williams and Wu . 


19931) . Therefore we decided to test the hypothesis whether the 


core of the Quadrantid meteor shower can be reproduced by 
meteoroid ejections around 1000 AD and 1500 AD. 

Fig. [TtI and [18] (upper panels) show the evolution of the 
descending nodes of meteoroids producing visual meteoroids, 
ejected from 2003 EHi in 1000 AD and 1500 AD respectively. 


1780-1882 



A-o, deg 


Figure 13l The simulated visual activity profile (blue boxes), superimposed 
over the observed activity profile (black dots) of the Quadrantid shower. Me¬ 
teoroids not approaching the Earth’s orbit within 0.01 AU are not plotted. The 
theoretical profile is a stack of 17 perihelion passages of 2003 EHi, correspond¬ 
ing to meteoroid ejections from 1780 - 1882 AD. The observed activity pro- 
file is an average act ivity profile of a few years between 1986 - 1992 (J2000) 
jRendtel et aliri993h . 
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Figure 14l The simulated radiant positions of the visual Quadrantids (black 
dots), along with the observed mean radiant positions and dispersions. The 
average radiant position is a stack of 17 perihelion passages of 2003 EHi, cor¬ 
responding to meteoroid ejections from 1780 - 1882 AD. Meteoroids which do 
not approach the Earth’s orbit within 0.01 AU are not plotted. The observed 
individual radiants and colors are the same as in Fig. ED 
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Figure 12: Evolution of the descending nodes of radar size meteoroids (upper panels), ejected from 2003 EHi in 1790 and 1796, until the present time. The red 
line marks the year of 1835 AD. The blue line marks the year of 2014, whereas the blue diamond con'esponds to the solar longitude of the peak of the radar shower. 
The lower panels show the simulated radiant positions (black dots). The green ellipse is centered at the mean radiant position as deduced from radar observations 
(CMOR) (see Table|4| Here the CMOR radiant is that deduce d from a 3D wavelet transform isolating the location of the peak average radiant between 2012 - 2014 
using the same technique as described in iBrown et al.N2010bh . In both figures, only simulated meteoroids approaching the Earth’s orbit within 0.01 AU are plotted. 
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Figure 15: The simulated radar activity profile (blue boxes), superimposed 
over the observed CMOR - measured activity profile (black dots) fro m single¬ 
statio n ech o observations following the same procedure described in iYe et all 
l2Q14h and jBrown and JonesL Il99^ . The theoretical profile is a stack of 15 
perihelion passages of 2003 EHi, corresponding to meteoroid ejections from 
1796 - 1886 AD. Only meteoroids approaching the Earth’s orbit within 0.01 AU 
are plotted. The observed activity profile is the average radar activity profile of 
the years 2002 to 2014. 


In this scenario we ejected 10"^ particles in a single perihelion 
passage of 2003 EHi in either 1000 AD and 1500 AD. 

Meteoroids ejected in 1000 AD and 1500 AD, show a bet¬ 
ter fit to the observed radiants positions (Fig. [T71 and [TSl lower 
panels) than ejection circa 1800 AD. However, in both cases 
the Quadrantid shower seems to first appear on the sky around 
1650 AD - too early to be consistent with first report of the 
shower around 1835. Furthermore, the same procedure was ap¬ 
plied to radar sized particles which also appeared as early as 
1650 AD. We also cannot uniquely distinguish which of the ob¬ 
served Quadrantids may be associated with the core and which 
are related to the activity wings which overlap in time. 
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Figure 16: The simulated radiants position of the radar Quadrantids (black 
dots), along with the observed (CMOR) mean radiant position and dispersion. 
The simulated radiants are from « 15 perihelion passages of 2003 EHi, cor¬ 
responding to meteoroid ejections from 1796 - 1886 AD. Only meteoroids ap¬ 
proaching the Earth’s orbit within 0.01 AU are plotted. 
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Figure 17l Evolution of the descending nodes of visual size meteoroids (up¬ 
per panels), ejected from 2003 EHi in 1000 AD until the present time. The red 
line mai'ks the year of 1835 AD. The blue line marks the year of 2014, whereas 
the blue diamond corresponds to the solar longitude (/!= 283.2°) of the peak of 
the radar shower. 

The lower panels show the simulated radiant position of meteoroids (black dots) 
superimposed over the observed mean radiant position, as measured by video 
and photographic techniques (see Fig. ED. In both figures, only simulated me¬ 
teoroids approaching the Earth’s orbit within 0.01 AU are plotted. 


7. Discussion and Conclusions 


We have used eight high-precision photographic Quadrantid 
orbits and integrated their orbits backward in time, along with 
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Figure 18: Evolution of the descending nodes of visual size meteoroids (up¬ 
per panels), ejected from 2003 EHi in 1500 AD, until the present time. The red 
line marks the year of 1835 AD. The blue line marks the year of 2014, whereas 
the blue diamond corresponds to the solar longitude of the peak of the radar 
shower. 

The lower panels show the simulated radiant positions (black dots) superim¬ 
posed over the observed mean radiant position, as measured by video and pho¬ 
tographic techniques (see Fig. IIU . In both figures, only simulated meteoroids 
approaching the Earth’s orbit within 0.01 AU are plotted. 


the assumed parent 2003 EHi, in order to constrain the most 
likely age of the core of the Quadrantid meteoroid stream. Out 
of the eight Quadrantids, five belong directly to the core of the 
stream having been observed within a day of the peak of the 
shower. In addition, five of the best quality radar Quadrantids 
detected by CMOR in 2013 were used as a complementary data 
set, in order to compare backward integrations of core radar¬ 
sized particles results with the visual counterpart of the stream. 

The most likely age of the core of the Quadrantid meteoroid 
stream, is a; 200 - 300 years, based on the backward integra¬ 
tions. This epoch was then used as a formation time window 
to test if the present observed characteristics of the radar and 
visual-sized particles in the ’core’ of the stream could be ex¬ 
plained by gas-drag ejection of meteoroids from 2003 EHi from 
this epoch. 

During the backward integrations, different values (0,10^^, 
10 “"^, 10“^) for p (solar radiation pressure to solar gravity) were 
considered in order to properly reflect the uncertainty in the 
mass of each individually observed Quadrantid. Our results did 
not show a noticeable difference as a function of different p 
values, which is not surprising given the timescale (200-300 
years) in our simulations. 

Our forward simulations indicate that an onset of cometary 
activity in 2003 EHi circa 1790 - 1796 provides a good match 
to all the observed constraints for the stream. In particular, the 
onset time of the shower in these forward simulations is near 
1830, consistent with historical data. Moreover, the location of 
the activity peak and the width of the core of the simulated ac¬ 
tivity profile are well reproduced in both photographic and radar 
simulations assuming ejection from this epoch. The broader 
(several day duration) weak activity surrounding the shower 
peak is not reproduced in our simulations, leading us to con¬ 
clude that while the core of the shower appears to be young, 
the wings are not. The simulated radiants of the radar and vi¬ 
sual showers also match observations well. We note that the 
simulated visual radiant is displaced approximately 1 degree 
from the mean observed radiant, but this is still reasonable as it 
lies within the broader s catter of the overall Quadrantid radiant 
( Jenniskens et al. . 1997h . 

In our forward modeling, for the m agnitude of the meteoroi d 
ejection speeds we used the model bv IBrown and JonesI jl998h . 
However, in order to check the robustness of our final results 
against different meteoroid ejection models, we repeated our 
forward simulations using the meteoroid ej ection model, re 


suiting from Crifo’s d istribu ted coma model (ICrifo et al.Lll995 


Crifo and Rodionovl 119970 . which generally yields slightly 


lower meteoroid ejection speeds compared to other models. In 
spite of the latter model being more sophisticated, the differ¬ 
ences in the final results were too small to change our conclu¬ 
sions. This perhaps should not be surprising given the time 
scales in our forward integrations (200-300 years). However, 
perhaps the difference in the final outcomes would be signif¬ 
icant for relatively young meteoroid trails (as young as a few 
orbital periods of the parent) which requires certain meteoroid 
ejection conditions to be fulfilled in order for the trail to inter¬ 
sect the Earth’s orbit. 

In a summary, our forward integrations are able to reproduce 
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the main features of the currently visible core portion of the 
Quadrantid stream using the formation epoch suggested from 
backward simulations («; 1800 AD years). The wider portions 
of the shower (lasting a few days) and the recently discovered 
(IBrown et al.L 1201 Obi) much longer radar activity (suggested as 
being a few months) are likely much older. 

At this time, it is difficult to disentangle the real na¬ 
ture of 2003 EHi, i.e a dormant comet or an asteroid. 
Could 2003 EHi have been active around 1700 - 1900 AD 
and gone undetected? It seems quite possible. Know¬ 
ing the absolute asteroidal magnitude of 2003 EHi (H=16.2, 
http://ssd.jpl.nasa.gov/sbdb.cgi) we calculated its 
visibility from Earth throughout the 19th century. Though 2003 
EHi reached almost 15th magnitude for a few brief periods in 
the 1800’s, it could easily have escaped detection, as its av¬ 
erage apparent magnitude was above 22. Of course, substan¬ 
tial cometary activity would have increased its brightness. Ac¬ 
cording to Kronk (Cometography volume 2, p ix, 2003), in the 
1800’s, small refractors (15-20 cm) were the most commonly 
used for comet observations. Larger ones were built towards 
the end of the century, with the 102 cm Yerkes telescope be¬ 
coming active in the late century, with ’’...the result was that 
early in the century comets were usually lost after having faded 
to magnitude 11 or 12, while at the end of the century comet s 
were being followed until near magnitude 16...” (Kronk, 20031). 
So it is plausible that even several magnitudes of brightening 
could have gone unnoticed at this time. Given the current inac¬ 
tive state of 2003 EHi it is entirely possible it was simply too 
faint, even if weakly active, to be detected as a comet. 

What might have triggered activity in a dormant parent? The 
orbit itself is stable in the century preceding our proposed start 
time, so it was not a dramatic orbital shift. Also, during the last 
few hundred years the perihelion distance of 2003 EHi has been 
consistently increasing, making it implausible that increasing 
solar heating triggered fresh activity. The ascending node of 
2003 EHi was near Jupiter in the 1700-1800’s, which makes 
collision with a main-belt asteroid unlikely, though a collision 
with a Jupiter Trojan remains a possibility. The descending 
node at the same time was near the Earth’s orbit, though no 
close approaches between any clones and Earth were recorded 
in our simulations. Thus a tidal encounter with Earth is not a 
likely cause either. 

We note that there was a relatively close approach between 
many of the clones of 2003 EHi and Jupiter in 1794, where the 
minimum distance reached 0.83 Hill radii. This was the closest 
approach to the giant planet since 1663 when it reached 0.69 
Hill radii. Thus, a reasonable explanation is that either tidal ef¬ 
fects from the encounter with Jupiter or a collision with a Trojan 
asteroid activated the parent for a time, with the activity subse¬ 
quently declining again to zero. However, the mechanisms of 
cometary activation and dormancy remain unclear and our pro¬ 
posed explanation is certainly not the only possibility. 
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